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ABSTRACT 

Submersibles  are  limited  in  depth  by  the  strength  of  the  pressure 
hull.  The  pressure  hull,  in  turn,  is  limited  in  strength  by  the  properties 
of  the  material  used  in  fabrication,  construction  techniques,  and  by  the 
veight  of  material  used  in  construction.  Once  the  ballast  weight  has  been 
taken  up  in  hull  veight,  a  limit  is  reached  in  depth.  To  go  deeper  "with  the 
same^payload,  the  hull  must  be  enlarged  or  size  may  be  held  constant,  the 
hull  strengthened  and  flotation  material  attached  to  support  the  extra 
veight.  Either  of  these  methods,  enlarging  the  hull  or  adding  flotation, 
increases  size  greatly  as  depth  is  increased. 

This  thesis  is  a  study  of  total  vehicle  size  versus  depth  for 
•weight -limited  deep  submersibles  vith  and  without  buoyancy  material 
attached. 


Thesis  Supervisor:  Philip  Mandel 

Title:  Professor  of  Naval  Architecture 
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I.  DISCUSSION 

A*  Volume-Limited  Submarines 

1.  Conventional  submarines  have  operated  routinely 
in  the  upper  layers  of  the  ocean  for  the  better  part 
of  a  century.  These  submarines,  designed  primarily 
as  military  vehicles,  are  large,  providing  the  neces- 
sary internal  volume  required  by  men  and  equipment. 
Weight  of  hull  and  equipment  is  less  than  the  sea 
vater  displaced,  thus  it  is  necessary  that  these 
vessels  carry  ballast  to  attain  neutral  buoyancy. 

2.  To  design  a  submarine  for  maximum  depth,  it  is 
necessary  that  the  hull  weight  be  increased.  This 
can  be  accomplished  with  a  volume-limited  submarine 
with  no  increase  in  size  by  trading  ballast  weight 
for  hull  weight.  The  depth  to  which  a  volume -limited 
submarine  can  be  designed  with  no  increase  in  size 

is  then,  a  function  of  the  ballast  it  carried.  Beyond 
that  depth  the  weight  of  hull  and  equipment  exceeds 
the  weight  of  sea  water  displaced,  and  the  vehicle 
becomes  weight -limited  requiring  additional  buoyant 
support  to  maintain  neutral  buoyancy. 
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B.  Deep  Submersibles 

1,   Oceano graph ers  require  small  laboratories  to.  carry 
them  to  the  depths  of  the  sea.  The  designer Ts  problem 
in  building  submarines  to  satisfy  the  needs  of  ocean- 
ography is  to  build  vehicles  capable  of  carrying  a 
scientifically  useful  payload  to  a  scientifically 
interesting  depth,  usually  the  maximum  depth  possible. 
Submarines  of  a  given  size,   -  limited  in  depth  by 
hull  veight  and  hull  materials,  have  several  avenues  of 
approach  by  which  their  depth  may  be  increased.  First, 
vehicles  may  be  made  larger  with  payload  fixed. 
Improved  pressure  hull  materials  may  be  used  or  flota- 
tion material  may  be  attached  to  the  vehicle  to  provide 
the  lift  necessary  to  restore  neutral  buoyancy.  The 
first  and  last  of  these  methods  result  in  large  craft  'which 
are  less  maneuverable,  more  complicated,  and  more  expensive 
to  build.  This  thesis  examines  and  compares  the  penalties 
paid  in  achieving  maximum  depth  by  the  methods  just  des- 
cribed utilizing  various  pressure  hull  materials. 
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II.  PROCEDURE 

A.  General  Approach 

1.  A  number  of  assumptions  are  made  concerning  the 
family  of  vehicles  to  he  considered,,  the  structural 
materials  available ,  the  desired  payload,  the  pro- 
pulsion characteristics  and  the  other  features  vhich 
vill  affect  the  buoyancy-weight  balance.  The 
assumptions  are  not  necessarily  optimum  for  the  range 
of  vehicles  considered ,  but  they  represent  an  average 
of  the  features  of  TRIESTE  II,  ALUMMAUT,  and  ALVBJ 
as  they  apply  to  the  family  of  vehicles  considered. 

2,  A  computer  program  was  developed  to  optimize  hull 
weight  for  a  given  depth.  The  program  is  described 
in  Appendix  A.  Using  the  computed  results  of  hull 
weight  versus  depth  as  a  base,  the  weight  of  payload, 
batteries,  dropable  ballast,  main  propulsion  units, 
and  "all  other"  variables  was  added.  This  resulted 
in  a  basic  vehicle  capable  of  proceeding  to  a  given 
depth  without  the  need  of  flotation  material.  At 
lesser  depths  additional  ballast  is  required  and  at 
greater  depths  added  lift  is  necessary. 
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3*     vehicles  of  various  size  vere  used  in  the  computation 
of  hull  weight  versus  depth.  Size  was  varied  in  two 
ways;  "by  holding  outside  diameter  constant  vhile  increas- 
ing length  and  "by  holding  the  length  to  diameter  ratio 
constant.  The  results  determine  changes  in  size  vith  an 
increase  in  depth  for  the  two  methods  of  increasing  dis- 
placement. 

k*     Finally ,  flotation  is  considered  as  a  means  of  hull 
support  at  increased  depth.  The  addition  of  flotation 
necessitates  an  addition  of  the  variables  mentioned 
above ,  "batteries,  "ballast,  etc.  This  reduces  the  effect 
of  the  flotation  material.  Curves  vere  constructed  for 
syntactic  foam,  probably  the  best  flotation  material  now 
available,  and  for  materials  of  lower  densities  which 
might  be  developed  in  the  not -too -distant  future.  These 
curves  take  into  account  the  weight  of  additional  equip- 
ment to  arrive  at  a  reduced  flotation  effectiveness. 
Entering  the  curves  vith  the  additional  lift  required  by 
a  basic  vehicle  at  a  given  depth,  one  determines  the 
volume  of  flotation  material  required,  weight  of  flotation, 
and  the  weight  of  the  added  hardware  necessitated  by  the 
additional  displacement.  The  result  of  primary  interest 
is  the  displacement  of  flotation  to  be  added  to  the  dis- 
placement of  the  basic  vehicle.  Total  displacement  is  then 
plotted  against  depth  for  comparison  of  size  with  similar 
curves  for  the  larger  pressure  hulls  considered. 
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B.  Assumptions 

1.  A  family  of  cylindrical  vehicles  with  hemispherical 
end  caps  is  considered. 

2.  A  "basic  vehicle  with  outside  diameter  of  seven  feet 
and  L/D  ratio  of  2.0,  where  L  is  the  length  including 
the  hemispherical  ends,  is  considered.  Hulls  with 
eight-foot  and  10-foot  diameters,  also  with  L/D  ratios 
of  2.0,  are  considered  for  comparison  as  are  vehicles 
of  seven -foot  diameter  with  L/D  ratios  of  2.5  and  3.0. 
Thus,  size  is  varied  in  two  ways;  "by  holding  L/D  constant 
and  "by  holding  D  constant. 

3.  Four  materials  are  considered  for  construction  as 
follows : 


Material 

Steel 

Steel 

Titanium 

Aluminum 


Yield  Strength 
80,000  p.s,i. 

150,000  p.s.i. 

120,000  p.s.i. 
60,000  p.s.i. 


Modulus  of  Elasticity 
30  x  10  p.s.i. 

30  x  10  p.s.i. 

6 
15  x  10  p.s.i. 


10  x  10  p.s.i. 


Density 
^90  lb. /ft.3 
^90  lh. /ft.3 
281  lb. /ft.3 
168  lh. /ft.3 


Glass 


Glass  was  also  included  for  comparison  in  the  computation 
of  hull  weights,  using  the  same  theory  applicable  to  metals. 
This  theory  may,  hut  probably  does  not  give  realistic  values 
for  the  weight  of  glass  hulls.  The  properties  used  for  glass 
are  as  follows; 

.  3Q0,000  p.s.i.       Ik  x  10  p.s.i.      184  lb. /ft.3 
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The  values  of  bull  weight  for  glass  are  tabulated 
in  Tables  I  -  III  but  are  not  plotted. 

4.  Power  requirements  are  based  on  a  sustained  speed 

of  five  knots  for  a  period  of  ten  hours.  The  auxiliary 
load  is  assumed  to  be  10  kilowatts. 

5.  Main  propulsion  power  is  assumed  proportional  to  the 
product  of  displacement  to  the  two -thirds  power  and  the 
speed  in  knots  cubed.  Since  the  speed  is  fixed  at  five 
knots,  power  is  proportional  to  displacement  to  the 
two -thirds  power. 

6.  Silver-zinc  batteries  with  the  following  properties 
are  assumed  for  the  power  source. 

Capacity  6.5   kilowatt  hours /cubic  foot 
23.O  pounds /kilowatt  hour^  ' 
7*  Payload  including  complement  and  instruments  is  held 

constant  at  2,000  lbs.  for  all  vehicles. 
8.  Minimum  volume  requirement  for  payload  and  equipment 
was  assumed  to  be  366  cu.ft.  which  is  the  internal 
volume  of  the  seven-foot  diameter,  L/D  =  2  pressure  hull. 


(l)  Keays,  K.,  "Parametric  Study  of  Two  Man  Search  and  Work  Submarines, 
General  Dynamics/Electric  Boat  contract  study,  page  24. 
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9.  Auxiliary  "buoyancy  is  assumed  initially  to  "be 
syntactic  foam  with,  a  density  of  k-k   lbs.  per  cubic 
foot.  Buoyant  materials  with  lesser  densities  of 
0.50  (32.1b.  per  cu.ft.)  and  0.30  (21.33  per  cu.ft. ) 
are  studied  for  comparison. 
10.  The  remaining  weights  which  do  not  fit  into  the 
categories  of  hull,  batteries,  or  payload  vary  with 
a  number  of  parameters.  Some  vary  with  displacement, 
some  with  depth,  some  with  speed  and  some  with  external 
or  internal  surface  area.  It  is  not  the  purpose  of 
this  paper  to  analyze  each  of  these  weight  groups,  but 
estimates  must  be  made  in  order  to  concentrate  on  the 
buoyancy  requirements.  It  is  assumed,  therefore,  that 
the  following  weight  categories  in  addition  to  hull, 
battery  and  payload  will  be  considered. 


Main  Propulsion  units  -  function  of  displacement  - 

3$  of  displacement 

Flotation  material  -  function  of  buoyancy  requirements 

and  depth 

Ballast  -  function  of  total  displacement  -  Qfo   of  displacement 

All  Other  -  for  the  purpose  of  this  thesis  a  function 
of  displacement  -  Qfo   of  displacement  C^). 


(2)  Wenk,  E.,  Jr.,  DeHart,  R.  C,  Mandel,  P.,  and  Kissinger,  R.,  Jr., 
"An  Oceanographic  Research  Submarine  of  Aluminum  for  Operation  to 
15,000  Feet",  Transactions  R.I.KJU,  Vol.  102,  No.  4.  Oct.  i960. 
Based  on  estimates  made  for  ALUMIKAUT,  page  15. 
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C.  Basic  Vehicle  Weights 

1.  A  computer  program  was  developed  "based  on  the  class 
notes  from  course  13.15;  M.  I.  T.,  "by  J.  Harvey  Evans. 
The  program  optimizes  frame  size  and  frame  spacing  to 
arrive  at  a  least  weight  solution  for  a  hull  of  a 
given  size,  constructed  of  a  given  material.  The  com- 
puter program  is  described  in  Appendix  A,  The  results 
are  tabulated  in  Tables  I  -  V. 

2.  Figure  VI  is  a  graph  of  battery  veight  versus  total 
displacement.  For  vehicles  above  30,000  lbs.,  battery 
veight  can  be  treated  as  a  constant  plus  a  percentage 
of  the  displacement  over  30,000  lbs.  with  an  error  less 
than  one -tenth  of  one  percent  of  total  hull  displacement 
at  all  times.  For  a  vehicle  to  ■which  flotation  material 
is  to  be  added,  the  additional  battery  weight  equals 
L.4  per  cent  of  the  additional  displacement. 

3«  Ballast,  propulsion  units,  and  "all  other"  also  are 
functions  of  displacement  and,  like  batteries,  these 
weights  may  be  treated  as  a  constant  which  is  part  of 
the  basic  vehicle,  plus  a  weight  equal  to  a  percentage 
of  the  displacement  of  the  buoyant  material  added.  Added 
weight  due  to  an  increase  in  displacement  amounts  to 
20.4$  of  the  added  displacement. 
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4.  Figures  II  -  V  are  plots  of  veight  versus  depth  and 
take  into  account  hull  veight,  batteries,  ballast, 
propulsion  units,  "all  other",  and  payload.  Weights 
plotted  in  Figures  II  -  V  are  exclusive  of  any  flotation 
material.  Displacement  is  also  plotted  in  Figures II  -  V 
for  the  differently  configured  vehicles.  Comparison  of 
displacement  to  veight  yields  the  ballast  or  the  added 
lift  required  for  neutral  buoyancy. 

D.  Buoyancy 

1.  The  buoyancy  of  a  vehicle  is  equal  to  the  veight  of 
salt  vater  displaced.  Flotation  material  must  be  added 
to  a  vehicle  vhich  veighs  more  than  it  displaces;  ballast 
must  be  loaded  aboard  a  vehicle  vhich  veighs  less  than 
its  displacement. 

2.  The  addition  of  flotation  material  necessitates  an 
addition  in  the  requirements  for  batteries,  dropable 
ballast,  main  propulsion  units,  and  "all  other".  Having 
linearized  battery  veight  as  a  function  of  added  dis- 
placement, all  of  the  above  items  may  be  treated  as  a 
single  veight  addition  proportional  to  the  flotation 
material  displacement. 

3.  The  lift  produced  by  flotation  material  is  also  a 
percentage  of  the  flotation  material  displacement.  The 
percentage  may  be  treated  as  the  effectiveness  of  the 
material.  Syntactic  foam,  at  kk   lbs.  per  cubic  foot 
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is  31* 35a  effective.  A  material  vith  a  density 
vith  respect  to  sea  vater  of  0.5  is  50$  effective 
and  one  vith  a  density  of  0.3  is  70$  effective. 

4.  The  effectiveness  of  the  flotation  material  Is 
reduced  by  the  additional  weight  described  in 
paragraph  2  of  this  section.  The  reduced  effective- 
ness is  the  flotation  material  effectiveness  less 
the  percentage  of  added  'weight.  For  the  assumptions 
made  in  this  study ,   the  reduced  effectiveness  of 
syntactic  foam  is  10. 9  (31.3  -  20.4);  of  the  0.5 
density  material  it  is  29.6;  and  of  the  0.3 
density  material  It  is  49.6.  Since  this  reduced 
effectiveness  is  the  true  indication  of  lift  contri- 
buted to  a  basic  vehicle,,  it  can  be  seen  that 
reducing  the  density  from  O.687  to  0*5  Increases  the 
effective  lift  almost  threefold.  Flotation  material 
displacement  is  reduced  to  one-third  of  the  initial 
requirement  and  added  batteries,  dropable  ballast , 
etc.,  is  only  one-third  as  much. 

5.  Figures  VTI  -  IX  are  graphs  of  flotation  material 
characteristics  and  veight  additions.  Entering  along 
the  ordinate  of  the  graph  vith  the  lift  required  by  a 
basic  vehicle  at  a  given  depth,  the  curves  provide  the 
following  information: 
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Curve  1  -  Lift  applied  to  the  "basic  vehicle 

Curve  2  -  Lift  applied  to  additional  "batteries, 
ballast,  etc. 

Curve  3  -  Total  lift  provided  by  the  flotation 
material 

Curve  k   -  Dry  -weight  of  flotation  material 
to  be  added 

Curve  5  -  Salt  water  displacement  of  flotation 
material  to  be  added 


These  curves  apply  to  any  of  the  vehicles  in  the  family 
of  vehicles  considered,  provided  that  flotation  material  is 
to  be  added  to  provide  additional  lift. 
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III.  PRESENTATION  OF  RESULTS 

Pressure  Hull  Sizes 
Figure  I 


\ 

! 


D  =  V ,   L/D  =  2.0,  2.5,  3.0 


D  =  8',  L/D  =  2.0 


D  =  10T,  L/D  =  2.0 


Displacement  of  Hulls 


Hull  Diameter 

L.O.A. 

L/D 

Displacement 
(Lbs.   Salt  Water) 

V 

lV 

17.5 

21.0 

2.0 

2.5 

3-0 

28,650 
37,300 
45,950 

8' 

16.0 

2.0 

42,900 

10' 

20.0 

2.0 

83,500 
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Table  I 
Computed  Cylinder  Weights  L/D  =1.5 


Collapse  Depth 


D 

3,4oo  ft. 

10,200  ft. 

17,000  ft. 

34,000  ft. 

Steel,  HY-80 

7' 

14,838  l"bs. 

23,800  lbs. 

45,400  lbs. 

8' 

8,3*5 

21,830 

34,900 

67,100 

10' 

15,600 

4l,600 

66,800 

127,400 

Steel,  HY-150 

7' 

8,600 

13,500 

25,200 

8' 

12,600 

19,700 

37,122   . 

10' 

23,900 

37,600 

71,000 

Titanium 

71 

6,120 

9,46o 

17,800 

8' 

3,620 

8,920 

13,900 

26,390 

10' 

6,760 

17,000 

26,600 

50,500 

Aluminum  60 

7r 

2,560 

10,600 

20,500 

8T 

3,74o 

9,780 

15,700 

30,300 

10' 

7,080 

18,600 

24,900 

58,400 

Glass 

7' 

1 

2,300 

3,400 

5,500 

8' 

3,500 

5,000 

8,100 

10' 

7,ioo 

10,000 

i5,4oo 
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Table  II 
Computed  Weight  of  2  Hemispherical  End  Caps 


Collapse  Depth 


D 

3,400  ft. 

10,200  ft. 

17,000  ft. 

34,000  ft. 

Steel,    HY-80 

7' 

2,300  lbs. 

7,300  lbs. 

11,900  lbs. 

22,700  lbs. 

8' 

3,700 

10,900 

17,800 

33,900 

10' 

7,200 

21,200 

34,700 

66,100 

Steel,   HY-150 

V 

1,320 

3,930 

6,500 

12,700 

8' 

1,980 

5,870 

9,700 

18,900 

10' 

3,860 

11,500 

18,900 

36,900 

Titanium 

7' 

1,050 

3,060 

5,000 

9,280 

81 

1,560 

4,550 

7,400 

13,900 

10' 

3,o4o 

8,900 

14,500 

27,100 

Aluminum  60 

V 

1,020 

3,030 

5,000 

9,700 

8' 

1,530 

4,500 

7,4oo 

i4,4oo 

10' 

3,000 

8,800 

14,500 

28,100 

. Glass 

7' 

250 

740 

1,23c 

2,400 

8' 

370 

1,100 

i,84o 

3,634 

10' 

730 

2,200 

3,600 

7,100 
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TaTole  III 
Computed  Pressure  Hull  Weights  L/D  =  2.5 


Collapse  Depth 


D 

3,4oo  ft. 

10,200  ft. 

17,000  ft. 

34,000  ft. 

Steel,  HY-80 

7T 

22,l4o 

35,700 

68,100 

8T 

11,850 

32,700 

52,700 

101,000 

10T 

22,800 

62,800 

101,500 

193,500 

Steel,   HY-150 

V 

12,530 

20,000 

37,900 

8J 

18,470 

29,400 

56,000 

10 T 

35,4oo 

56,500 

107,900 

Titanium 

7T 

9,380 

i4,46o 

27,080 

8J 

5,180 

13,^70 

21,300 

40,250 

10T 

9,800 

25,900 

4i,ioo 

77,6oo 

Aluminum  60 

7T 

3,58o 

15,600 

30,200 

8' 

5,270 

14,300 

23,100 

44,700 

10' 

10,080 

27,400 

44,400 

86,500 

Glass 

7r 

3,o4o 

4,630 

7,900 

8' 

4,600 

6,84o 

11,730 

10 T 

9,300 

13,600 

22,500 
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Table  IV 
Computed  Pressure  Hull  Weights  L/D  =  2.0 


Collapse  Depth 


D 

3,400  ft. 

10,200  ft. 

17,000  ft. 

34,000  ft. 

Steel.  HY-80 

7r 

17,200  lbs. 

27,800  lbs. 

52,300  lbs. 

8T 

9,ioo 

25,400 

41,100 

78,700 

10 ' 

17,600 

49,000 

79,200 

151,100 

Steel,  HY-150 

7r 

9,700 

15,500 

29,500 

8' 

14,300 

22,800 

43,700 

10' 

27,400 

44,000 

84,200 

Titanium 

7r 

7,150 

11,300 

21,200 

8' 

4,000 

10,500 

16,700 

31,400 

10? 

7,500 

19,700 

32,200 

60,500 

Aluminum  60 

V 

2,730 

12,000 

23,400 

8' 

4,100 

11,000 

18,100 

34,600 

10' 

7,700 

21,200 

34,4oo 

67,000 

Y 
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Table  V 


Computed  Pressure  Hull  Weights  l/D  =3.0 


Collapse  Depth 

D 

3,400  ft. 

10,200  ft. 

17,000  ft. 

34,000  ft. 

Steel,   HY-80 


7r 

27,100  lbs. 

43,700  lbs. 

81,900  lbs. 

8' 

14,500 

40,000 

64,400 

123,500 

10' 

28,000 

76,800 

123,700 

236,100 

Steel,   HY-150 


7' 

15,400 

24,500 

46,300 

8' 

22,700 

35,900 

68,500 

10' 

43,300 

69,100 

131,500 

Titanium 


1  " 

71 

11,200 

17,600 

33,ooo 

8T 

6,500 

16,600 

26,300 

49,500 

10 T 

12,200 

32,000 

50,500 

95,500 

Aluminum 


7' 

4,400 

19,100 

37AOO 

8' 

6,700 

17,500 

28,800 

54,8oo 

10' 

12,400 

33,600 

54,300 

105,900 
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Table  VI 

Weights  Added  to  Bare  Hulls  to  Arrive 
at  Basic  Vehicle  Weights 


(Figures  II-V) 

Diameter 

7' 

V 

7' 

l/d 

2.0 

2.5 

3.0 

Battery  Wt. 
(Figure  VI ) 

3,300 

3,550 

3,750 

Ballast,  Etc. 
(19$  x  Displ. ) 

5,^50 

7,100 

8,7^0 

Payload 

2,000 

2,000 

2,000 

Total 

io,750 

12,650 

1M50 

8* 

10* 

L/D 

2.0 

2.0 

Battery  Wt. 

3,6oo 

4,300 

Ballast,  Etc. 

8,150 

15,900 

Payload 

2,000 

2,000 

Total 

13,750 

22,200 
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Figure  XII 
Sample  Computer  Solution 


Material  -  Aluminum  60 
Collapse  Depth  -  If, 000  feet 
Hull  Diameter  -  7T 


Shell 


1.46' 


T 

8.68' 


11.7" 


L—  9.36'L-J 


3.02" 


-25.2" 


5.02" 

_J_ 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Figure  X  displays  the  results  required  "by  section  I  of  this  thesis. 
Overall  vehicle  size  is  plotted  against  depth  for  hulls  constructed  of 
four  different  materials.  The  effect  of  increasing  depth  by  changing 
hull  size  or  "by  adding  flotation  material  vith  pressure  hull  size  con- 
stant is  compared  vith  the  following  results.  Increasing  hull  size  by 
varying  vehicle  length  provides  greater  depth  than  building  a  larger 
vehicle  vith  L/D  constant  for  HY-80  and  Titanium.  For  aluminum  the 
opposite  is  true  and  for  RY-150  both  methods  have  equal  effect.  Therefore,, 
the  results  vith  respect  to  a  preferred  method  of  increasing  size  are 
inconclusive.,  one  method  having  nearly  the  same  effect  as  the  other. 
However,  the  addition  of  flotation  material  is  to  be  preferred  over 
changing  vehicle  size  by  either  method  for  all  but  HY-80,  assuming 
syntactic  foam  is  the  material  used.  For  depth  increases  greater  than 
about  5.> 000  feet,  a  smaller  vehicle  results  if  flotation  material  is 
added.  And,  of  course,  the  cost  of  a  given  volume  of  flotation  material 
vill  be  much  less  than  the  cost  of  increasing  hull  size  by  an  equal 
amount. 

Recommendation:  Deep  submergence  pressure  hulls  should  be 
built  large  enough  to  satisfy  minimum  volume  requirements.  If  added 
lift  is  necessary,  it  should  be  provided  by  flotation  material. 
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3.  The  effect  of  improving  the  characteristics  of  flotation  materials 
is  graphically  demonstrated  in  Figure  X.  If  a  flotation  material  could 
he  produced  with  a  density  of  0.50,  there  -would  he  no  question  ahout  the 
preferred  method  of  increasing  depth. 

Recommendation:  As  flotation  material  of  lover  density  is 
developed,,  its  use  should  he  considered  in  the  design  of  deep  submersibles. 
Costs  will  he  cut  and  smaller,  less  complex  deep  suhmersibles  will  result. 

C.  Although  the  results  concerning  a  preferred  method  of  changing  vehicle 
size  were  inconclusive,  it  can  he  said  that  there  is  probahly  an  optimum 
hull  configuration.  With  the  computer  program  developed,  one  could  arrive 
at  hull  proportions  for  a  least  weight  vehicle  with  a  given  displacement. 
At  first  glance  it  might  appear  that  this  would  always  result  in  a  sphere, 
hut  with  conventional  framing  or  sandwich  hull  construction,  this  probably 
will  not  be  the  case.  The  program  can  also  be  used  to  determine  optimum 
frame  configuration,  and  it  appears  to  be  a  very  promising  tool  in  evalu- 
ating current  design  practices  as  they  apply  to  vehicles  designed  for 
greater  depths. 

Recommendation:  The  computer  program  described  in  Appendix  A 
should  be  used  to  continue  the  study  of  hull  configuration  started  in  this 
paper  and  it  could  be  used  as  well  for  a  study  of  the  parameters  of  sub- 
marine design  and  their  interaction  on  one  another. 

D.  Figure  XI  is  a  graphic  illustration  of  the  results  obtained  by  using 
the  methods  outlined  in  the  notes  from  course  13*15  aacl  from  earlier  cal- 
culations based  on  unstiffened  cylinder  construction.  Some  difference 
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exists  bet-ween  the  weight  of  unstiffened  cylinders  and  the  weight  of 
the  stiffened  cylinder  as  predicted  "by  the  program  in  Appendix  A, 
It  is  natural  that  a  difference  exist,  but  in  some  cases  the  unstiff- 
ened  cylinder  is  limited  to  a  depth    38$  less  than  that  of  a  stiffened 
cylinder  of  equal  weight.  This  large  difference  is  difficult  to  explain 
and  it  may  he  an  indication  that  the  notes  on  which  the  computer  program 
is  based  may  give  incorrect  collapse  depths  when  considering  deep  sub- 
mersibles,  or,  that  the  computer  program  contains  some  flaw.  However, 
the  program  was  tested  using  as  inputs  data  from  relatively  shallow- 
depth  submersibles,  and  it  duplicated  almost  exactly  the  calculations 
which  resulted  when  those  vessels  were  designed. 

Recommendation':  The  computer  program,  probably  the  most 
valuable  by-product  of  this  thesis,  should  first  be  checked  once  again 
to  ensure  that  the  formulae  have  been  translated  into  computer  language 
properly,  and  then  the  output  should  be  compared  with  experimental  results 
to  determine; 

1.  the  validity  of  the  theory  with  respect  to 
conventional  steel  construction, 

2.  the  applicability  to  other  metals,  and 

3.  the  application  of  the  theory  to  vessels 
with  t/D  ratios  approaching  0.1. 

Having  checked  the  program  and  its  general  applicability,  it  will  be 
possible  to  determine  the  necessity  of  modifying  the  basic  theory  for 
metals  other  than  steel  or,  more  probably,  for  deep  submergence  vehicles 
with  very  thick  shells. 
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APPEITOIX  A 

Computer  Design  Program  to  Optimize  Submarine  Hull  Weight,  "based  on 
class  notes  of  Course  13.15,  M.  I.  T.,  by  J.  Harvey  Evans,  Professor 
of  Naval  Architecture. 

1 .  Background 

a.  Failure  of  a  stiffened  shell  occurs  in  one  of  three 
modes.  The  first  occurs  if  the  shell  "buckles  "between 
the  frames  while  the  frames  retain  their  circularity. 
This  is  referred  to  as  elastic  or  plastic  instability. 
The  second  mode  of  failure  involves  yielding  of  the 
shell.  In  this  mode  the  shell  remains  circular  "but 
the  yield  strength  of  the  material  is  exceeded.  The 
third  mode  of  failure  occurs  when  the  shell  and  frames 
fail  together  and  this  is  known  as  general  instability. 
Ideally,  for  optimizing  strength  with  least  weight,  a 
shell  should  fail  in  all  three  modes  simultaneously. 
This  program  optimizes  three  parameters,  shell  thickness, 
frame  area,  and  frame  spacing  in  an  attempt  to  find  the 
least  weight  solution. 
.  b.  In  such  a  solution  it  would  be  ideal  to  optimize  the 
dimensions  of  the  frame  as  well  as  the  other  parameters 
mentioned,  but  this  would  have  introduced  another  degree 
of  complexity  to  an  already  complex  analysis.  Instead, 
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a  frame  typical  of  those  presently  used  in  submarine 
construction  "was  selected  as  a  parent  and  all  frames 
are  scaled  from  the  parent  member.  The  cross  section 
of  the  parent  member  is  shown  in  Figure  A-l  below. 

All  dimensions  are     with  respect  to  Tt  r.  web 

•^        w  ' 

thickness. 


Shell 
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c.  Having  selected  a  parent  frame  member,,  it  is  only 
necessary  to  settle  on  the  length  of  cylinder  between 
bulkheads  and  the  mean  cylinder  diameter  in  order  to 
applythe  formulae  presented  in  the  set  of  notes  refer- 
enced above.  The  general  method  of  approach  outlined 


-  35  - 


is  that  of  selecting  a  set  of  scantlings  and  checking 
their  applicability.  An  understanding  of  the  paper  is 
necessary  for  a  full  understanding  of  the  program  and 
its  capabilities.  However,  lack  of  such  understanding 
will  not  prevent  successful  use  of  the  program, 
d.  The  program  is  based  on  the  solution  of  six  equations. 
These  are: 

1.  Hoop  stress  formula  which  predicts  a  pressure 
at  which  shell  yield  occurs. 


Per  =  2(t/D)  cry 


This  formula  is  used  to  determine  shell  thickness 
for  a  contemplated  design.  Shell  thickness, 

Per  x  D 

Experiments  indicate  that  the  predicted  thickness 
may  be  reduced  in  some  cases  by  as  much  as  20$. 
The  computer  program,  however,  completes  two 
solutions,  at  predicted  't1  and  at  95$  of  'tr  but 
a  minor  modification  to  the  program  allows  'tr  to 
be  reduced  to  any  given  percentage  of  the  value 
predicted  by  the  hoop  stress  formula,  or  it  allows 
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'tT  to  be  reduced  in  increments  until  a  least 
weight  solution  is  found  which,  satisfies  all 
other  criteria  hut  ignores  the  hoop  stress 
prediction, 
2*  Formula  92a  modified  predicts  plastic  insta- 
bility failure. 
2t 

-         m  J i 


1  +  :i0-85  -  B    1  +  0.l48(\2  -  ef 
1  +  3 


The  critical  pressure  it  predicts  must  he  equal 
to  or  greater  than  collapse  depth  pressure,   (see 
page  42  for  a  list  of  symbols  and  their  meanings . ) 
3.  WenkTs  formula  predicts  collapse  pressure  based 
on  the  circumferential  stress  at  the  shell  mid- 
thickness  midway  between  frames. 

2ta 


P   -  -^  x 


cr     D 


i^hl 


4.  Lunchick  predicts  "yield  pressure"  based  on  a 
three  hinge  failure  mechanism. 


2tc 
P  =~-£  x 

y   d 


Vx2  -  xz  +  z2 


This  formula  is  normally  within  3  -  5$  of  the  collapse 
pressure  predicted  by  Wenk. 
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5.  Bryant  has  taken  a  collapse  prediction  by 
Kendrick  and  put  it  into  a  useful  form.  This 
formula  is  "based  on  the  number  of  complete  cir- 
cumferential waves  formed  by  the  plating  and 
frames  buckling  as  a  unit.  It  is  solved,,  normally 
for  2,    3,   or  h  waves  and  the  first  critical  mode 
of  failure  (i.e.,  the  one  for  which  the  critical 
load  is  smallest)  is  the  value  sought. 


cr 


2Et 


k 


x 


m 


,2   m    _  v ,  2    2  N2 
(n  +-jj--l)(n  +  m  ) 


+  8<n  -  ^  EL 


Formulae  2,   3,   H-  and  three  variations  of  5,   predict 
critical  pressures.  It  is  necessary  to  solve  these 
formulae  and  then  select  that  answer  which  is  lowest 
as  the  critical  one  in  the  case  under  consideration, 
6.  The  last  formula  used  predicts  the  stress  in  the 
frame  flange  due  to  the  combined  effects  of  total 
load  including  shear  and  bending  due  to  out  of 
roundness.  This  formula  is  based  on  the  critical 
pressures  predicted  by  Bryant's  formula j  thus,  it 
is  also  solved  three  times  in  the  program. 


max 


PR 


1  + 


0.8501 
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1  +  3 
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P   -  P, 
cr    1 
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The  highest  values  predicted  by  this  formula 
must  not  exceed  the  yield  strength  of  the  material. 
Thus,  the  highest  value  is  selected  "by  the  program 
as  the  critical  value. 

e.  The  program  solves  for  shell  thickness  using  the  hoop 
stress  formula  and  then  computes  six  solutions  for  critical 
pressure.  The  lowest  of  the  six  is  termed  the  critical 
pressure  for  solution  of  the  problem.  Then  three  solutions 
of  maximum  stress  are  made,  and  the  highest  of  these  is 
selected  as  the  critical  stress.  At  this  point  the  program 
computes  the  amount  hy  which  critical  pressure  must  he 
changed  to  agree  with  collapse  pressure,  and  the  amount  hy 
"which  the  highest  stress  must  be  altered  so  It  does  not 
exceed  yield  strength  of  the  material.  If  the  program  can 
he  made  to  alter  frame  size  and  frame  spacing  until  these 
values  agree;  critical  pressure  with  collapse  pressure,  and 
maximum  stress  with  yield  stress,  then  for  a  given  shell 
thickness,  a  least  weight  solution  will  he  possible. 

f .  The  solution  is  found  by  solving  the  problem  with  frame 
area  incremented  a  -unit  amount  and  with  frame  spacing  incre- 
mented a  unit  amount.  These  two  solutions  yield  changes  in 
the  critical  values  of  pressure  and  stress;  that  Is,  they 
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give  the  partial  differentials  of  pressure  and  stress 
■with,  changes  in  frame  area  and  frame  spacing.  It  is 
then  a  simple  matter  to  assume  that  linearization  is 
possible  and  solve  two  equations  in  two  unknowns  for 
the  number  of  unit  changes  to  "be  made  in  frame  area,  A, 
and  frame  spacing,  SP. 

|~dA  +  ~|dSP  =  (P  __     -P  ..) 
dA      oSP      v  collapse    crit' 


||dA  +  J|cLSP=(S..--S   ..) 

oA      oSP       v  yield    crrir 


Solving  for  the  two  unknowns,  dA  and  dSP,  the  original 
estimate  of  frame  area  and  frame  spacing  can  be  changed. 
The  problem  is  solved  again  with  the  new  values  of 
frame  area  and  frame  spacing  and  if  the  critical  pressure 
is  within  a  given  percentage  (the  program  specifies  kfo 
on  the  safe  side)  and  if  the  critical  stress  is  also  safe 
and  within  a  set  tolerance  (these  tolerances  can  be 
changed  if  a  higher  degree  of  accuracy  is  desired),  the 
program  will  compute  weight  of  compartment  including  shell 
and  frames.   (One  word  of  warning.  In  order  to  ensure  a 
satisfactory  solution  in  all  cases,  it  was  necessary  to 
make  frame  spacing  continuous.  Therefore,  the  number  of 
frames,  equal  to  compartment  length  over  frame  spacing 
may  not  be  a  whole  number.  It  will  be  necessary  to  select 
the  next  higher  number  of  frames  to  ensure  a  safe  solution 
with  a  slightly  higher  weight  than  that  computed. ) 


g.  If ,  at  the  end  of  the  first  iteration,  the  values  of 
predicted,  collapse  pressure  and  predicted  stress  are  not 
vithin  tolerance _,  the  program  will  continue  to  cycle 
through  the  series  of  calculations  improving  the  solution 
until  it  is  vithin  tolerance.  The  program  then  increments 
shell  thickness  to  95$  of  its  predicted  value  and  solves 
the  problem  again.  The  weight  determined  is  compared  with 
the  weight  determined  with  the  original  value  of  lt r  and 
the  "least  weight"  is  printed  "before  the  problem  goes  on 
to  the  next  set  of  data. 


-  in 


LIST  OF  SYMBOLS 

Some  of  the  symbols  in  the  following  list  refer  to  variables 
■which  are  used  to  simplify  the  basic  formulae.  When  such  variables  do 
not  have  explicit  definitions _,  reference  is  made  to  a  formula  number, 
figure  number ,    or  page  number,  vhere  the  symbol  is  used  in  the  notes 
by  Professor  Evans. 

A  Cross  sectional  frame  area 

A,  Area  of  one  frame  plus  shell  under  the  web 

b  Web  thickness  of  the  ring  frame 

B  Page  18 

c  Formula  31 

D  Outside  diameter  of  hull 

D      Diameter  of  hull  to  mid-shell 
m 

E  Young's  Modulus 

F  Page  31 

H  Figure  8 

I  Moment  of  inertia  of  a  T  section 

L,  Frame  Spacing 


m 


n 


3.1^-  D  -  2L  (rounded  off  to  the  nearest  whole  number) 

2,    3,  or  h  -   no.  of  lobes  for  the' three  modes  of 
failure  by  general  instability 


P      Predicted  collapse  pressure 


P      Pressure  at  which  yield  occurs 
cr 
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P   .,   Lowest  collapse  pressure  of  the  6  pressures  calculated 

P  Yield  pressure  -  same  as  P 
V  cr 

P  Pressure  at  collapse  depth 

R  Radius  of  hull  to  mid-shell 
m 

S    ,   Highest  stress  of  the  three  computed 

SP  Frame  spacing 

t  Shell  thickness 

x  Page  3^ 

z  Page  3^ 

P  Page  18 

e  Page  9 

A  Slenderness  ratio  -  Figure  k 

YS  Yield  stress 
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2 .   Input 


a.   Input  consists  of  the  length  and  mean  diameter  of  the 
cylinder  under  consideration,  the  characteristics  of  the 
material  to  "be  used  in  construction,  including  yield 
strength,  Young's  modulus  and  density  of  the  material, 
the  collapse  depth,  allocable  out -of -roundness  of  the 
hull,  and  the  "best  estimate  available  of  frame  area  and 
frame  spacing. 

Definition 
Area  of  frame 
Collapse  depth 
Compartment  length 
Mean  diameter  of  hull 
Allowable  out -of -roundness 
Density  of  material 
Young's  modulus 
Frame  spacing 
Yield  strength 
(Delta,  allowable  out-of -roundness  of  the  hull,  should 
normally  be  limited  to  one -half  inch  or  half  the  thick- 
ness of  the  shell.  For  heavy  shell  plating  this  may  be 
increased,  but  for  shell  plating  less  than  one-inch 
thick  Delta  will  be  altered  by  the  program  to  one-half 
the  thickness  of  the  shell  plating. ) 


Inputs 

Units 

A 

2 
in 

CD 

feet 

CL 

inches 

D 

inches 

Delta 

inches 

Dens 

lb/in3 

E 

p.s.i. 

SP 

inches 

SY 

p.s.i. 
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b.  All  input  for  a  single  problem  is  placed  on  one  card 
in  the  following  form: 


1  A    7   SY     17   E 


27  D   33  CD  39  CL   5Dens  51SD  55Delta 


The  numbers  indicate  the  card  space  number  at  which  each  item 
should  start.  Sample  input  follows: 

70 60000. 30000000.-360.—  4oo.-54o.  —  .2875-30.  -.50- 

The  sample  input  indicates  the  following: 


Frame  area 

Yield  Strength 

Young's  Modulus 

Mean  diameter  of  hull 

Collapse  depth 

Compartment  length 

Density  of  material 

Frame  spacing 

Allowable  out-of -roundness 


70  sq..in. 
60,000  p.s.i. 
30  x  10  p.s.i. 
360  inches 
4,000  feet 
540  inches 
.2875  lb.  per  in.3 
30  inches 
0.50  inches 
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An  input  deck  ready  for  the  computer  vill  consist  of  the 
following  cards : 

Meaning    Space   Number 

7 

M3999-350^FMS,DEBUG,   1,   1,    500,    75         Peterson 

XEQ 

BINARY     (or  FORTRAN) 

Peterson     (this  card  included  only  if  a 
fortran  deck  is  used) 


1 

Identity 

■K 

Execute 

*• 

Binary 

* 

Comment 

c 

First  binary  or  fortran  card  of  program 


Last  card  of  program 

END 

DATA 

First  data  card 
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MINIMUM  HULL  WEIGHT  PROGRAM 


10 


20 
30 
40 
50 


15 


60 

64 

65 
66 

67 


DIMENSION  PCAM(6)  ,PCW(6)>PCL(6),PCB1(6)  ,PCb2(6)  ,PC33(6 
1) ,  SMK1 (6) »SMK2 (6 ) ,SM.<3(6 )  ,  P  (  1 0 ) , S ( 4 )  , W ( 40 )  » PD I FS ( 1 0 ) 
2,SDIFS( 10) 

READ  1340,A,SY,E,D,CD,CL»DEi\S,SP,DT.LTA 

PCD=CD*64./144. 

PRINT  1010, SY 

PRINT  1040, CD 

PRINT  1090, PCD 

T=CD*64.*D*0.5/ (144.*SY) 

IF(T-0.25)20,30,30 

T=0.25 

IF(T-1.0)40,50,50 

DELTA=T/2. 

CONTINUE 

SY=0.98*SY 

CD-1.02*CD 

PRINT  1020 


PRINT 
PRINT 
PRINT 
PRINT 


tE 
,D 
,CL 
,DENS 


CALCULATIONS 


1030 

1050 

1060 

1080, DELTA 
PCD=1.02*PCD 
LR  =  1 
L=l 
LO=l 
LP  =  1 
NP=1 
NQ=1 

PRELIMINARY 
FN=CL/SP+1. 
IF(A)  64,64,65 
PRINT  1390, A 
GO  TO  10 
IF(SP)66,66,67 
PRINT  1400,  SP 
GO  TO  10 
LQ  =  1 

TW=SQRTF( A/19. 18) 

THETA=18.2*SP/ (D*SQRTF( 100.*T/D) ) 
SINH1=( 2.7182  8**THETA-1./2.7182  8**THETA) /2.0 
C0SH1=( 2.71828**THETA+1./2.7182  8**THETA: /2.0 
GK= ( S INH1-SI NF ( THETA) )/ (SI NH1+SINF(THETA) ) 
GN=(COSHl-COSF(THETA) ) / ( Si NH1+S I NF ( THET A ) ) 
THA2=THETA/2. 

SINH2=( 2.7182  8**THA2-1./2.71828**THA2 >/2. 
COSH 2=  (  2.7  182  8**THA2+1./2.71828**THA2  )  /2. 
GH=-(3.*SINH2*COSF(THA2)+COSH2*SINF(THA2) ) / ( S I NH1+S I NF 
KTHETA)  ) 

Q=-( SINH2*COSF(THA2 ) +COSH2*S I NF ( THA2 ) ) / ( SI NH1+SINF ( THE 
1TA)  ) 
GT=-(SIISH2*COSF(THA2)-COSH2*SINF(THA2)  )/(SINHl+SINF(TH 
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1ETA)  ) 

BDA=( ( SP/ D ) **2 /(T/D)**3)**0. 2 5*SQRTF( SY/E) 

EPSL=(0.45*SY*D)/(E*T) 

B=TW*T/(  A+TW*T  ) 

BETA=( ll.*GN*T**2 )/(SQRTF( 100.*T/D )* ( A+TW*T ) ) 

BM=3.1416*D/  (2.-N-CL) 

IPM=XFIXF( BM+0.5 ) 

BM=FLOATF( IBM) 

F=1.+2.*Q*{ 0.85-B) /( l.+BETAJ 

X=l .+GH*( 0.8  5-B )/ ( 1,+BETA) 

Z=0.5+3.6  2*GT*{0.8  5-B)/ { 1.0+BETA) 

EL*1.57*SQRTF(D*T/2.  ) 

IF (EL-SP)80»80,70 
70    EL  =  SP 

8  0    ERTIA1=102.2*TW#*4*(3.*EL*T/(EL#T+19.18*TW**2)*(0.125* 
1T/TW+1.43 )**2+0.79  ) 

AB=A+TW*T 

C= C43.5*TW**3+ (T*EL)*( 8 . 0*TW+0 . 5*T ) )/ ( 19.18*TW**2+EL*T ) 

PC AM ( L)=2.*T*SY/(D*( l.+GH* ( (0.8  5-B)/ ( 1 . +BETA )  ) ) * ( 1 .+0  . 
1148*(  (8DA**2-EPSL)**2 )  )  ) 

NO  =  l 

P(NO)=PCAM(L) 

NO=NO+l 

PCW(L)=2.*T*SY/(D*SQRTF(F**2-F/2.+0.2  5) ) 

P(NO)=PCW( L) 

NO=NO+l 

PCL(L)=2#0*T*SY/(D*SQRTF(X**2-X*Z+Z*Z) ) 

P(NO)=PC^(L) 

NO=NO+l 

BN  =  2. 

PCB1(L)=2»*E*T*BM**4/(D* (BN**2+( .5*BM)**2-1« )*(BN#*2+B 
1M**2)**2 )+8.*(BN**2-l. )*E*ERTIA1/ (D**3*SP) 

P(N0)=PCB1 (L) 

NO=NO+l 

BN  =  3. 

PCB2 (L )=2.*E*T*BM**4/(D*lBN**2+( .5*BM)**2-1. )*(BN**2+d 
1M**2  )**2)+8.*(BN**2-l.  )*E*ERTIA1/  (D**3*SP.) 

P(NO)=PCB2( L) 

NO=NO+l 

BN  =  4. 

PCB3( L)=2.*E*T*BM#*4/{D* (BN#*2+< «5*BM)**2-1« )*(BN**2+d 
1M**2)**2)+8.*(BN**2-1,)*E*FRTIA1/ (D**3*SP) 

P(NO) =PCB3 ( L) 

BN  =  2. 

SMK1 ( L)= (PCD*(0.5*D+0,5*T)*TW/ (AB*( 1.0+BETA) ) }*( 1.0+(0 
1.85*BETA)/B) 
2+E*C*(BN**2-1.0)*DELTA*PCD/ ((0.5*D)*-*2*(PCB1(L)-PCD) ) 

NO  =  l 

S(N0)=SMK1 (L) 

NO=NO+l 

BN  =  3. 

SMK2 ( L)=( PCD*( 0.5*D+0.5*T)*TW/ ( AB* ( 1 . 0+BETA ) ) ) * ( 1 . 0+ ( 0 
1.85*BETA)/B) 
2  +  E*C*(BN**2-l«0-)*DELTA*PCD/  ((0.5*D)**2*(PCB1  (L)-PCD)  ) 

S(NO)=SMK2(L) 
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N0=N0+1 

BN  =  4. 

SMK3 ( L ) = (PCD*( 0.5*D+0.5*T)*TW/ (AB*(1.0+BETA) ) )* ( 1 .0+(0 
1.85*BETA) /B) 

2+E*C*(BN**2-1.0)*DELTA*PCD/  (  ( 0.5*D )**2*( PCSl ( L)  -PCD)  ) 
S(N0)=SMK3(L) 
:      COMPUTE  VOLUME  OF  MATERIAL  IN  ONE  "OMPARTMENT  LENGTH 

V=0.7  8  54*CL*(  ( D+T/2. )**2-( D-T/2. ) *~2 ) +0 . 78  5  4* ( (D-T/2.  ) 
l**2-(D-T/2.- 

2  5.93*TW)**2  )#TW*FN+0.7854*FN*6.4*TW*(  (  D-T/2  .-5.  93#TW  )  * 
3*2-(D-T/2.-8.0*TW)**2) 
:      COMPUTE  WEIGHT  OF  MATERIAL  IN  ONE  COMPARTMENT  LENGTH 
W( L) =V*DENS 
L  =  L  +  1 

GO  TO  (90,165) ,LP 
90    DO  110  N=l,6 

PDIFS(N)=PCD-P(N) 
110   CONTINUE 

PDIFCR=PDIFS(1 ) 
DO  120  N  =  2,6 

IF(PDIFS(N)-PDIFCR) 120*100,100 
100   PDIFCR=PDIFS(N ) 

NP  =  N 
120   CONTINUE 

122   IF(0.02*PCD-ABSF(PDIFCR) )130,125,125 
125   LQ=LO+l 
130   DO  140  N=l,3 

SDIFS(N)=SY-S(N) 
140   CONTINUE 

SDIFCR=SDIFS( 1 ) 
DO  160  N=2,3 

IF (SD I FS(N)-SDIFCR) 150,160,160 
150   SDIFCR=SDIFS(N) 

NQ  =  N 
160   CONTINUE 

162  IF( 0.02*SY-ABSF(SDIFCR) )164, 163,163 

163  LO=LQ+l 

GO  TO  ( 164,164,200) ,LQ 

164  LP=2 
PCRIT=PCj-PDIFCR 
SCRIT=SY-SDIFCR 

165  GO  TO  (  170,180,190)  ,LO 
170   AT=A 

A=0.96*A 

LO  =  2 

GO  TO  60 
180   DELPE1=PCRIT-P(NP) 

DELS1=SCRIT-S(NQ) 

A  =  AT 

SPT=SP 

SP=1.04*SP 

LO  =  3 

GO  TO  60 
190   DELPE2=PCRIT-P (NP) 

DELS2=SCRIT-S(NQ) 
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200 


SP 
DO 
DO 
PR 
PR 
A  = 
SP 
GO 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PP 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 
PR 


=  SPT 
T=DELS 
P=SDIF 
Y=DOP/ 
X=(PDI 
(l.+O. 
=(1.-. 
TO  15 


2-DELPE2*DELSl/DELPEl 

CR-PDIFCR*DELS1/DELPE1 

DOT 

FCR-PRY*DELPE2)/DELPEi 

04*PRX )*A 

04-PRY  )*SP 


INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 
INT 


10 
10 
11 
11 
U 
11 
12 
12 
12 
12 
12 
11 
11 
U 
U 
11 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
45 
11 
11 
11 
11 
11 
11 
11 
U 


00 
70 
40 
"0 
oO 
70 
00 
10 
20 
30 
40 
00 
10 
20 
80 
90 
0 


A 

FN 
PC 
PC 
PC 
PC 
PC 
PC 
SM 
SM 
SM 
SP 
T 

TW 
V 
W( 
,THE 

♦  SIN 

♦  COS 
»GK 
,GN 
,THA 

♦  SIN 

♦  COS 

♦  GH 

♦  Q 

♦  GT 
>BDA 

♦  EPS 

♦  5 

♦  BET 
tBM 

♦  F 

♦  X 

♦  Z 

♦  EL 

♦  ERT 

♦  AB 

♦  C 
30,  P 
30,  SC 
30, PD 
30.  S 
JO, 
30, 
30, 
30, 


AM(L-1 ) 

W (L-l ) 
L (L-l  ) 
Bl (L-l  ) 
B2(L-1  ) 
B3(L-1  ) 
Kl (L-l  ) 
K2(L-1  ) 
K3(L-1 ) 


L-l  ) 
TA 

HI 
HI 


2 
H2 

H2 


IA1 


CRIT 

RIT 

IFCR 

DIFCR 

ELPE1 

ELPE2 

ELS1 

ELS2 
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PRINT  1130,  DOT 

PRINT  1130,  OOP 

PRINT  1130,  PRY 

PRINT  1130,  PRX 

GO  TO  ( 210,220,250) ,LR 
210   WLEAST=W(L-1 ) 

PRINT  1370,  WLEAST 

LR  =  2 

TR  =  T 

T=0.9  5*T 

GO  TO  15 
220   IF(WLEAST-W(L-1  )  ) 240 , 240 »2 30 
230   WLEAST=W(L-1 ) 
240   PRINT  1370, WLEAST 

T=1.04*TR 

LR  =  3 

GO  TO  15 
250   IF(WLEAST-W(L-1 ) )270, 270*260 
260   WLEAST=W(L-1 ) 

PRINT  1370,  WLEAST 

T=1.04*T 

GO  TO  15 
270   PRINT  1370, WLEAST 

GO  TO  10 
450   FORMAT( FL0.5 ) 

1000  FORMAT( lHltlOHFRAME  AREA > 1 9X , 1H= , FlO . 2 , 3X , 7HSQ.  IN.) 
1010  FORMAT( 1H1,14HYIELD  STRENGTH , 15X , 1H= , Fl 0 . 2 , 3X , 6HP . S. I . ) 
1020  FORMAT ( 1X,14HY0UNGS  MODULUS , 12X , 1H= , F13 . 2 , 3X , 6HP . S . I . ) 
1030  FORMAT( 1X,21HMEAN  DIAMETER  OF  HULL  ,  8X  ,  1H= , Fl 0 . 2 , 3X , 6H I 

NCHES) 
1040  FORMAT(  IX  ,  1 4HCOLLAPSE  DEPTH  ,  1  5X  ,  1H=« »  F  10  .  2  ,  3X  ,  4HFEET  ) 
1050  FORMAT ( IX , 16HBULKHEAD  SPAC I NG , 13X , 1H= ,F1 0 . 2 , 3X , 6H I NCHE 

S) 
1070  FORMAT ( IX, 30HNO.  FRAMES  BETWEEN  BULKHEADS  =,F10.2,3X,6 

HFRAMES) 
1060  FORMAT ( IX , 1 9HDENS I T Y  OF  MATER  I AL , 10X , 1H= , Fl 0 . 4 , 3X , 16HL 

3.  PER  CU.  I 
INCH) 
1080  FORMAT( IX , 30HALLOWABLE  OUT-OF-ROUNDNESS    =,F10.2,3X,6 

HINCHES) 
1090  FORMATC IX , 30HPRESSURE  AT  COLLAPSE  DEPTH    =,F10.2,3X,6 

HP.S.  I  .  ) 
1100  FORMAT { 1X,13HFRAME  SPAC I NG , 16X , 1H= , Fl 0 . 2 , 3X , 6H I NCHES ) 
1110  FORMAT( 1X,15HSHELL  TH I CKNESS , 14X , 1H= , Fl 0 . 2 , 3X , 6HI NCHES ) 
1120  FORMAT* IX, 13HWE6  TH I CKNESS , 16X , 1H= , FlO . 2 , 3X , 6HI NCHES ) 
1130  FORMAT(F20.5) 
1140  FORMAT( 1HO,30HCOLLAPSE7PRESS.  FROM  92AMOD.  =,F10.2,3X, 

6HP.S.  I.  ) 
1150  FORMAT( IX , 30HCOLLAPSE  PRESS.  FROM  WENK     =,F10.2,3X,6 

HP.S. I . ) 
1160  FORMAT( IX ♦ 30HCOLLAPSE  PRESS  FROM  LUNCH  I CK  =,F10.2,3X,6 

HP.S. I  .  ) 
1170  FORMAT( 1X.30HCOLL  PRESS  FR  BRYANT,  2  LOBES= , F 10 . 2 , 3X , 6 

HP.S.  I.  ) 
1180  FORMAT( 1X,30HVOLUME  OF  ONE  COMPT  MATERIAL  =,F10.2,3X,1 
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2HCUBIC  INCH) 
1ES) 
1190  FORMAT  (  1X.21HWEIGHT  OF  COMPARTMENT  »  8X  »1H=»F10.  2  >.'  X»6HP 

OUNDS ) 
1200  FORMAT ( IX t30HCOLL  PRESS  FR  BRYANT,  3  LOBES=  ,  F 10  .  2 , 3X , 6 

HP.S. I . ) 
1210  FORMATt 1X.30HCOLL  PRESS  FR  BRYANT,  4  L03ES= , F 10  .  2  ,  3X , 6 

HP.S. I . ) 
1220  FORMATt IX, 30HMAX  STR  FR  BUSHIPS,   2  L03ES  =,F10.2,3X,6 

HP.S. I . ) 
1230  FORMATt IX, 30HMAX  STR  FR  b'JSHlPS,   3  LOBES  =,F10.2,3X,6 

HP.S. I . ) 
1240  FORMAT ( IX, 30HMAX  STR  FR  BU'HIPS,   4  LOBES  =,F10.2,3X,6 

HP.S. I  .  ) 
1340  FORMAT ( F6.2»2F10.0,3F6.0*F6.4,F4.0,F4.2) 
1370  FORMATt 1X.30HCOMPARTMENT  OF  LEAST  WEIGHT   =,F10.2,3X,6 

HPOUNDS) 
1390  FORMATt IX, 30HACTION  TERMINATED,  NEGATIVE  A= , F10 . 2 , 3X , 7 

HSQ.  IN.) 
1400  FORMAT(1X,30HACTION  TERMINATED,  NEGAT.  SP  =,F10.2,3X,6 

HINCHES) 
END 


NOTE 

FORMAT  STATEMENTS  REQUIRE  A  CARD  FORM  72  CHARACTERS 
IN  WIDTH.  IN  ORDER  TO  FACILITATE  A  NARROWER  PRINTOUT 
THOSE  FORMAT  STATEMENTS  LONGER  THAN  60  CHARACTERS  WERE 
BROKEN  AND  THE  EXTRA  CHARACTERS  PRINTED  ON  THE  FOLLOWING 
LINE  WITH  AN  EXAGERATED  INDENTATION.  A  USER  COPYING  THE 
PROGRAM  SHOULD  TYPE,  FOR  EXAMPLE,  LINE  1400  AND  THE  NEXT 
INDENTED  LINE  ON  THE  SAME  CARD. 
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COMPUTER  PROGRAM  SYMBOLS 


A         Cross-sectional  area  of  a  frame 

AB        Cross-sectional  area  of  a  frame  plus  the  area  of 
shell  under  the  "web 

AT  Same  as  A 

B  Page  18 

EDA  Lambda  -  Figure  k 

BETA  Page  18 

BM  m  (see  formula  symbols) 

BN  n  (see  formula  symbols) 

C  Formula  31 

CD  Collapse  depth 

CL  Compartment  Length 

C0SH1  cosh  e 

C0SH2      COSH  6/2 

D         Hull  diameter 

DELPE1     P   . ,  -  P  ,,  (after  A  is  incremented) 
crit    crit  x  ' 

DELPE2     P   , .  -  P  ..  (after  SP  is  incremented) 
crit    crit  v  ' 

DELS1      S   . .  -  S   .,  (after  A  is  incremented) 
crit    crit  x  ' 

DELS2      S  . .  -  S  ..  (after  A  is  incremented) 

DELTA     Allocable  out -of -roundness  of  the  hull 

DENS      Density  of  construction  material 

DOT,  DOP   Intermediate  values  in  the  solution  of  two 
differential  equations 

E         Young rs  modulus 
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Hit 

EPSL 

ERTIA1 

F 

FN 

GH 

GX 

GN 

GT 

P 

PCAM 
PCB1 

PCE2 

PCB3 

PCD 

PCL 

PCRIT 

PC¥ 

PDIFCR 

PDIFS 

PRX 

PRY 

Q 


Effective  length  of  shell  plating 

e,   page  9 

Koment  of  inertia  of  a  T  section 

Page  31 

No.  of  frames  in  one  compartment  length 

Figure  8 

Figure  8 

Figure  8 

Page  3^ 

Predicted  collapse  pressure  =  PCAM,  PCW,  PCL,  PCB1,  PCB2, 
PCB3 

Predicted  collapse  pressure  from  modified  formula  $2k 

Predicted  collapse  pressure  from  Bryant rs  (Kendrick) 
formula  (1st  failure  mode) 

(2nd  failure  mode) 

(3rd  failure  mode) 

Pressure  at  collapse  depth 

Predicted  collapse  pressure  from  Lunchick 

Lowest  collapse  pressure  of  the  6  predicted 

Predicted  collapse  pressure  from  Wenk 

Highest  of  the  6  PDIFS  values 

=  PCD  -  P  (six  of  these  calculated) 

No.  of  increments  by  -which  SP  must  be  changed  in 
reaching  an  optimum  solution 

No.  of  increments  by  vhich  A  must  be  changed  in  reaching 
an  optimum  solution 

Page  31 
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S  Predicted  stress  =  SMKL,  SxMK2,  SMK3 

SCRIT  Highest  stress  of  the  3  predicted 

SLIFCR  Lovest  of  the  3  SDIFS  values 

SDIFS  SY  -  S  (three  calculated) 

SINH1  SUffi  e 

SINH2  SINH  0/2 

SMK1  Predicted  stress  (2  lobe  mode  of  failure) 

SMK2  Predicted  stress  (3  lobe  mode  of  failure) 

SMK3  Predicted  stress  (h   lobe  mode  of  failure 

SP  Frame  spacing 

SPT  Same  as  SP 

SY  Yield  stress 

T  Shell  Thickness 

THA2  0/2 

THETA  0  (Figure  8) 

TW  Web  thickness 

V  Volume  of  material  in  shell  and  frames  over  one 
compartment  length 

W  Weight  of  material  in  shell  and  frames  over  one 
compartment  length 

WLEAST  Weight  of  the  lightest  compartment  of  those  considered 

X  Page  3^ 

Z  Page  3*f 
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i|.  Output 

a.  Given  the  following  inputs  for  a  section  of  pressure 
hull: 

1.)  Diameter 

2.)  Length  "between  transverse  "bulkheads 

3. )  Collapse  depth 

K.)  Construction  material 

5.)  Tolerance  to  which  the  construction 

material  will  he  fabricated 

the  computer  program  produces  recommended  values  for; 
1.)  Shell  thickness 
2.)  Frame  size,  and 
3.)  Frame  spacing. 

Wrihthe  above  information,  it  is  possible  to  compute 
displacement  of  the  cylinder  and  weight  of  material  used 
in  construction  of  the  hull.  The  computer  program  does, 
in  fact,  calculate  weight  of  hull  since  it  was  decided 
that  "least  weight"  would  be  the  criterion  for  optimum 
design. 

b.  But  the  process  of  design  involves  more  than  one  answer 
to  one  problem.  Normally,  a  designer  will  study  several 
vehicles  with  various  characteristics  and  he  will  attempt 
to  gather  as  much  information  as  possible  about  each  one. 
It  is  not  enough  to  know  that  a  computer  program  has 
recommended  certain  values  for  a  certain  vehicle.   Each 
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of  the  recommended  values  must  be  checked  for  accuracy } 
logic,  and  applicability  to  the  planned  use  and  the 
construction  methods  currently  available.  There  is 
always  the  possibility  that  a  computer  program  ■will 
recommend  as  optimum  a  set  of  values  impossible  to  adhere 
to  in  actual  construction. 

To  aid  a  designer  in  comparing  vehicles  and  to  aid  in 
analysing  the  ccuiputer  p^o^ram  itself,  or  possibly  to 
help  in  studying  the  design  parameters  of  submarine  con- 
struction, the  computer  program  listed  in  this  section 
prints  out  a  large  number  of  values  computed  in  the 
preliminary  calculations  as  well  as  the  final  answers 
to  the  main  problem.  Once  a  user  has  gained  familiarity 
with  the  program  he  "will  probably  eliminate  much  of  the 
printed  output  by  'withdrawing  from  the  input  deck  those 
PRINT  statements  referring  to  the  unnecessary  values. 
He  may  also  add  some  cards  or  some  additional  calculations 
■which  will  help  him  in  his  design  or  analysis. 

Referring  to  the  sample  output  (page  oO  )_,  the  following 
values  are  printed  out  after  each  solution  at  a , given 
thickness. 
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Computer  Output 
Page  One 

1.  Yield  strength  of  the  material 

2.  Collapse  depth 

3.  Pressure  at  collapse  depth 
h.  Young's  modulus 

5.  Hull  diameter 

6.  Compartment  length 

7.  Density  of  the  material 

8.  Allocable  out -of -roundness  of  the  hull 

Page  Two 

1.  Prame  Area 

2.  Number  of  frames 

3-8.  Critical  pressure  predicted  by  formula  92a- 
modified,  Wenk,  Lunchick,  and  three  modes 
by  Bryant 

9-11.  Three  values  of  yield  strength  predicted  by  the 
formula  in  section  l-c-6  of  this  appendix. 

12.  Frame  spacing 

13.  Shell  thickness 

14.  Web  thickness 

15.  Volume  of  material  in  the  shell  and  frames 

16.  Weight  of  material  in  the  shell  and  frames 

17-38.  Preliminary  calculation  values  beginning  vith  "THETA" 

and  ending  ■with  "C"  as  listed  on  the  fourth  page  of  the 
computer  program. 

39-50.  Values  used  in  the  solution  of  the  two  differential 
equations  when  determining  the  required  changes. 

51.   "Least  ■weight"  of  those  compartments  considered. 
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d.  Since  thickness  is  first  decreased  to  95$  of  the 
predicted  value  and  then  increased  to  104$  of  that 
predicted  by  the  hoop  stress  formula,  the  values  on 
page  two  will  he  printed  at  least  three  times.   If, 
as  may  happen  in  some  cases,  the  solution  with 
thickness  increased  produces  the  "least  weight" 
solution,  thickness  will  continue  to  he  increased  by 
increments  of  4$  until  a  "least  weight"  solution  is 
reached  and  passed.   (Each  of  these  incremented 
solutions  will  also  require  one  page  of  printed 
output.)  This  series  of  events  would  occur  if,  by 
increasing  thickness,  frame  size  and  the  number  of 
frames  could  he  decreased  enough  to  achieve  an 
overall  reduction  in  weight.  After  a  certain  number 
of  thickness  increases,  however,  weight  will  begin 
to  increase  at  which  time  the  program  will  terminate. 
Experience  indicates,  however,  that  the  "least 
weight"  solution  will  occur  when  "T"  is  reduced  to 
95$  of  the  thickness  predicted  hy  the  hoop  stress 
formula,  and  if  the  hoop  stress  prediction  is  ignored, 
"least  weight"  will  occur  at  about  70$  of  the 
predicted  thickness. 
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YI5LC  STRENGTH 

CCLLAPSE  DEPTH 

PRESSURE  AT  CCLLAPSE  DEPTH 

YCUNGS  MODULUS 

MEAN  DIAMETER  CF  HULL 
8LLKHEAD  SPACING 
DENSITY  CF  MATERIAL 
ALLOWABLE  CUT-CF-RGUNONE SS 
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P.S.I . 

1530. 

FEET 

68C.:: 

P.S.I. 

3  3  0  0  C  C  • 

P.S.  I. 

214.70 

INCHES 

416.00 

INCHES 

.2  840 

LB.  PER 

.50 

INCHES 

CU.  INCH 
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FRAME  AREA 
N  C .  FRAMES 


BETWEEN  BULKHEADS  = 


8.26 
19.56 


SQ.  IN. 
FRAMES 


CCLLA 

CCLLA 

CCLLA 

CCLL 

CCLL 

CCLL 

MAX  S 

MAX  S 

MAX  S 

FRAME 

SHELL 

WEB  T 

VCLUM 

WEIGH 


PSE7 
PSE 
PSE 
PRES 
PRES 
PRES 
TR  F 
TR  F 
TR  F 
SPA 
THI 
HICK 
E  CF 
T  CF 


PRES 
DRES 
PRES 
S  FR 


FR 

FR 
6U 
BU 
BU 


CING 

CKNE 

NESS 

ONE 

COM 

2 

5 

6 

1 
1 
1 
1 


18 

131 

9 

4 

663 

55986 

-6 

347 

1 

7 

-219 

-87 

1104 

3239 

2 

6 

COMPARTMENT  C 


S  .  f 

S .  FROM 
S  FROM 
BRYANT 
BRYANT 
BRYANT 
SHIPS, 
SHIPS, 
SHIPS, 

SS 

COMPT 

PARTMEN 
.49059 
.99276 
.07  562 
.81633 
.041 27 
.24  530 
.59  3 "5 
.88091 
.5:168 
.34727 
.  19287 
.65292 
.  14018 
.  G  9  0  2  8 
.53025 

.ccoco 

.85262 
.89355 
.64815 
.  18325 
.  12429 
.03498 
.70594 
.14374 
. 13652 
.42158 
.75732 
.  3  8  3  6  C 
.50674 
.65967 
.05908 
.70824 
.C3876 
.93203 
. 10457 
F  LEAST 


9? 

W  E 
LUN 
,  2 
,  3 
i  4 

2 

3 

4 


A  '•'  0  0.   : 
CHICK   : 

LOSES 

L  0  B  E  S : 

LOBES' 

LOBES  '■ 

LOBES  ■■ 

LOBES  : 


MATERIAL  = 
T 


70C 

853 

791 

47?3 

15  58 

2207 

41390 

46662 

54042 

22 

1 

228702 
64951 


'2 
14 
74 
31 
19 
23 
5C 
0  6 
24 
41 
25 
66 
30 
45 


S. 
S. 


P.S.  I 
P.S.I 
P.S 

s 
s 

s 
s 


INCHES 
INCHES 
INCHES 

CUBIC  INCHES 
POUNDS 


'EIGHT 


63257.29    POUNOS 
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